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Wound healing is a complex process in which 
a tissue's individual cells have to be orches- 
trated in an efficient and robust way. We in- 
tegrated multiplex protein analysis, immunohistochemical 
analysis, and whole-slide imaging into a novel medium- 
throughput platform for quantitatively capturing prolif- 
eration, differentiation, and migration in large numbers 
of organotypic skin cultures comprising epidermis and 
dermis. Using fluorescent time-lag staining, we were able 
to infer source and final destination of keratinocytes in 
the healing epidermis. This resulted in a novel extending 



shield reepithelialization mechanism, which we confirmed 
by computational multicellular modeling and perturba- 
tion of tongue extension. This work provides a consis- 
tent experimental and theoretical model for epidermal 
wound closure in 3D, negating the previously proposed 
concepts of epidermal tongue extension and highlight- 
ing the so far underestimated role of the surrounding 
tissue. Based on our findings, epidermal wound closure 
is a process in which cell behavior is orchestrated by a 
higher level of tissue control that 2D monolayer assays 
are not able to capture. 



Introduction 

In human skin wound healing, reepithelialization is the most es- 
sential part, as the tissue's primary objective is to quickly rees- 
tablish barrier function (Martin, 1997; Singer and Clark, 1999; 
Friedl and Gilmour, 2009). The individual cells of the skin are 
orchestrated to behave in such a way that skin integrity is rees- 
tablished in an evolutionarily proven, most robust way (Singer 
and Clark, 1999). It is highly challenging to design experiments 
capturing how this orchestration actually takes place. Although 2D 
monolayer experiments are ideal for analyzing individual cellular 
functions such as migration mechanistically on the single cell 
level, wound healing cannot be reduced merely to cell migration 
(Farooqui and Fenteany, 2005; Soderholm and Heald, 2005; 
Liang et al., 2007). Thus, for understanding wound healing, the 
analysis of the orchestration of the individual processes taking 
part in wound healing has to be performed. This can only be un- 
dertaken in 3D wound-healing models, which have to be system- 
atically and quantitatively characterized. The goal is hereby to 
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derive consistent computational models helping to uncover high- 
level tissue functions as well as to understand the roles of indi- 
vidual cellular processes in tissue repair. In the sense of Noble 
(2006), it is the question of how a repair function at the higher bi- 
ological scale of the tissue is actually realized by the lower scale 
of the single cell level. 

Choosing this systems biological approach can be ex- 
pected to provide answers to several open questions of wound 
closure. A central question debated in the literature in skin 
wound healing is, for example, the mechanism of the creation 
and extension of the epidermal tongue. Two reepithelialization 
mechanisms were postulated so far. The first is the leap-frog or 
rolling mechanism in which migrating suprabasal cells roll over 
leading basal cells and dedifferentiate to form new leaders 
(Krawczyk, 1971; Paladini et al., 1996). The tractor-tread or 
sliding mechanism postulates that layered keratinocytes move 
forward in a block (Radice, 1980; Woodley, 1996). A variant is 
the model of Usui et al. (2005) in which suprabasal cells migrate 
out of the wound, thereby outnumbering the basal cells. It has 
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Figure 1 . The excisional 3D in vitro wound model using precast EFT cultures and DEs. (A) Hematoxylin and eosin-stained paraffin sections of the EFT 
cultures. The epidermis (£) displayed basal, spinous, granular, and cornified keratinocyte layers (C), and the dermis (D) distributed fibroblasts. Broken lines 
denote dermal-epidermal junction. (B) Precast EFT skin cultures were wounded twice with a 2-mm biopsy punch. (C) After wounding, cultures were placed 
into a 1 2-well ThinCert insert. (D) Cultures were halved with one punch on every tissue sample and cut until the wound was microscopically visible. Sections 
were marked as layer 1 (0 pm). Based on that, further defined sections in layer 2 (150 urn) and in layer 3 (300 pm) were produced. (E) Histology of the 
wound cultures. Hematoxylin and eosin-stained sections showed punched wounds with clearly defined wound margins (arrowheads). The underlying 
DE was seamlessly attached, providing the wound matrix (dotted line shows basal membrane). (F) Hematoxylin and eosin reepithelialization kinetics of 
the wound cultures over a period of 1 0 d. Broken lines denote wound margins. (G) The relative wound closure was quantified as the ratio of the migration 
distance (length of the EET; /eet) to the size of the wound in the respective section (/wound). The data shown are from two independent experiments; n = 8 
(day 0 [d0]), 10 (day 1), 8 (day 2), 10 (day 3), 9 (day 4), 4 (day 5), and 3 (day 6); 100 slides were analyzed in total; data are means ± SEM; **, P < 
0.01 ; ***, P < 0.001, Student's t test. Bars, 500 pm. 



up till now been unclear whether one of these mechanisms is 
correct and how such a mechanism is functionally embedded in 
the environment of the wound. The latter issue points to the 
question of the contributions of the intact surrounding tissue, 
which has been largely neglected so far and thus warrants a sys- 
tematic analysis. Both aspects, tongue extension and the intact 
tissue of the wound, are linked to and realized by tightly regu- 
lated spatiotemporal processes of proliferation, migration, and 
differentiation, finally leading to reestablishment of the intact 
epidermal 3D morphology of the skin (Gurtner et al., 2008). 

To build a consistent mechanistic model of wound closure, 
we set up a dedicated technical analysis pipeline comprising 3D 
organotypic wound models, standardized immunohistology, 
fluorescent whole-slide imaging, image analysis, multiplex pro- 
tein analytics, and computational systems biological modeling. 
We applied our pipeline on large numbers (92) of 3D organo- 
typic full-thickness skin wound models comprising keratinocytes 
and fibroblasts, which we tracked in time by a novel two-step 
time-lag fluorescence staining. This allowed us to dissect the 
epidermal 3D wound-healing process spatiotemporally in cell 
proliferation, migration, and differentiation and to derive the 



extending shield mechanism (ESM), a consistent theory of how 
these three processes are intertwined leading to the incremental 
and robust closure of human wounds. 

Results 

The organotypic skin wound model shows 
a concentric closure with rates similar to 
the in vivo situation 

To construct a reproducible wound model, commercially obtain- 
able epidermal full-thickness (EFT) cultures, containing epidermis 
and dermis, were punched and adhered to self -fabricated dermal 
equivalents (DEs; Fig. 1 , A and B). With a handling time of 5 min per 
3D culture involving wounding, attaching to the DE, and plac- 
ing the wound culture into inserts, our experimental setup enables 
a considerable medium throughput (Fig. 1 C). For subsequent 
histological sectioning, we developed a defined tissue-sectioning 
method (see Materials and methods) to ensure consistent and 
comparable wound diameters for analysis (Fig. 1 D). Similarly 
as intact epidermis, in vivo, the unwounded epidermis of the cul- 
tures contained about five cell layers, in which basal keratinocytes 
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showed a typical palisade structure, characterized by cubical cell 
shape and continual flattening with progressive differentiation, 
finally building up the stratum corneum. All EFT cultures seam- 
lessly attached to the DE, providing an in vivo-like matrix for 
reepithelialization surrounded by clearly defined wound margins 
(Fig. 1 E). 1 d after wounding, keratinocytes had started migration 
(Fig. 1 F). At this point, the extending epidermal tongue (EET) was 
visible in the form of two cell layers at the margins and one cell 
layer at the leading edge. 4-5 d after wounding, keratinocytes had 
covered the whole wound area, forming a neo-epidermis, which 
displayed a thickness of about three cell layers near the wound 
margins and one cell layer at the center of the wound. From 
day 4-10, the neo-epidermis thickened and formed a multilayered 
epithelium. 1 d after wounding, 36% of the former wound area 
was closed (Fig. 1 G). At day 4 after wounding, reepithelialization 
of 100% of the wounded area was achieved, equaling complete 
wound closure, leading to a reepithelialization rate of 338 um/d 
being in the physiological range (Laplante et al., 2001). 

Keratinocyte-fibroblast cross talk ensures 
an in viva-like wound-healing environment 

Keratinocyte-fibroblast cross talk and the underlying paracrine 
growth factor signaling are essential to coordinate wound 
healing-associated processes such as proliferation and migration. 
To show that our wound model shares important components 
of the epidermal-dermal cross talk with native skin wounds, 
we measured selected growth factors driving the reepithelial- 
ization mechanism. We collected supernatants of wounded and 
unwounded organotypic cultures and analyzed IL-la, IL-3, 
IL-6, IL-8, MCP-1, and granulocyte colony-stimulating fac- 
tor (G-CSF) concentrations before (day 0), as well as 1 and 3 d 
after wounding (Fig. 2, A and B). Most abundant were the proin- 
flammatory cytokines IL-6, IL-8, and MCP-1, showing a mean 
concentration of >20,000 pg/ml. G-CSF (> 12,000 pg/ml) was 
also present in high concentrations followed by IL-3 (>800 pg/ml) 
and IL-la (>500 pg/ml). Although IL-la, IL-3, IL-6, and G-CSF 
peaked 1 d after wounding, IL-8 showed an continuous increase. 
Elevated levels of MCP-1 concentration were only found 3 d after 
wounding, assuming a delayed tissue reaction. Based on the 
cytokine measurements, the underlying paracrine signaling 
network was drawn in Fig. 2 C. 

Keratinocyte proliferation occurs in a 
concentric traveling wave after the EET 

To determine the spatiotemporal profile of keratinocyte pro- 
liferation during reepithelialization, we divided histological 
sections of the wound cultures in 10 regions equaling five con- 
centric rings and analyzed their epidermal proliferation index 
(EPI) as the fraction of proliferating basal cells over 10 d 
(Fig. 3 A). Comprehensive computational image analysis of the 
digital slides yielded an EPI baseline of 13% after averaging 
the EPI at day 0 and 10. Emanating from this baseline, the whole 
culture reacted with a proliferative burst at day 1 with an EPI 
of <50% within the wound and the surrounding tissue (Fig. 3 B). 
At this time point, the wound itself still showed only few in- 
migrated keratinocytes. With progressing time and accompa- 
nying increasing migration, the proliferative activity declined 



depending on the distance to the wound margin. 3 d after wound- 
ing, the EPI significantly decreased in the surrounding tissue, 
while remaining continuously high with 42% in the wound. The 
EPI then normalized until day 6 in the surrounding tissue as well 
as in the wound area. Based on phosphorylated histone 3 stain- 
ing, the mitotic index was calculated showing a spatiotemporal 
profile equal to the EPI (Fig. SI, A-C). We performed system- 
atic morphometric analysis of the unwounded tissue regions 
with image processing of histological sections. Right after 
wounding, we observed only a limited increase in the thickness 
of the CK10 - cell layers despite nonthickening of CK10 + layers 
(Fig. SI D). From day 4 on, the height of the CK10~ layer was 
even below those of day 0. In light of the massive and sudden 
proliferation response of the whole tissue, we conclude that the 
newly produced cells must have been used to establish a mas- 
sive concentric cell stream toward the wound. Collectively, the 
results show an immediate burst in proliferation in a concentric 
pattern around the punch wound, which traveled as a wave to- 
gether with former and newly produced cells into the direction 
of the punch wound. 

Collective polarization of basal 
keratinocytes in the intact tissue 
toward the wound 

We analyzed the surrounding tissue for signs of collective po- 
larization, indicating a potential cell stream toward the wound, 
indicating in turn collective migration. Therefore, we analyzed 
nuclear shape, positioning, and deformation as well as Golgi 
complex translocation during wound closure. We separated the 
organotypic skin cultures into six regions of equal size accord- 
ing to the wound distance and quantified the major axis (see 
Materials and methods) of every nucleus of the basal cells by 
computational image processing (Fig. 4 A). Unwounded cultures 
showed a distribution of nuclear orientations with a tendency to- 
ward 90°, in which 0° represents a horizontal orientation. After 
12 h, we observed in all regions (R1-R6) of wounded cultures 
a distribution that significantly differed from rotation angles of 
nuclei in unwounded skin (Fig. 4 B). Thereby, the rotation was 
strongest in the inner rings R3 and R4 and quantitatively di- 
minished with increasing distance from the punch. 36 h after 
wounding, this trend progressed to an even stronger rotation in 
which many cells took a fully horizontal orientation with their 
major axis parallel to the basal membrane. Only a neglectable 
(0.6%) fraction of cells still showed a vertical (90°) orientation. 
This observed shift of nuclear rotation from vertically to hori- 
zontally orientated nuclei was in line with the translocation of 
the Golgi complex, a cellular polarization marker, which shifted 
from a position underneath the nuclei in unwounded cultures 
to a position behind the nuclei in an opposed direction to mi- 
gration in wound cultures (Fig. S2). We further qualitatively 
measured the position of the nucleus relative to the cell bound- 
aries stained with E-cadherin (Fig. 4 C). Cell nuclei showed 
a center position as expected in unwounded cultures, whereas 
in wounded cultures, nuclei moved rearward, opposed to the 
distantly located punch wound. At the same time, nuclei of ro- 
tated keratinocytes showed an elongated, cigarlike morphology. 
Nuclei of keratinocytes in unwounded cultures showed a round 
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Figure 2. Paracrine signaling of keratino- 
cytes and fibroblasts in the excisional 3D 
in vitro wound-healing model. (A and B) 

Mean concentrations of growth factors mea- 
sured in the supernatant of the 3D in vitro 
wound-healing model. The growth factors are 
divided into two groups, based on the maxi- 
mum mean concentration: >1 0,000 pg/ml (A); 
<1,000 pg/ml (B). (C) The paracrine signal- 
ing of keratinocytes and fibroblasts under 
wound-healing conditions is shown. Arrows 
mark dependencies between the cell types. 
Dashed lines represent possible dependen- 
cies, which could not be measured in the 3D 
in vitro wound-healing model caused by the 
lack of leukocytes. Growth factor-secreting 
cell types and the impact of the factors on cell 
behavior are listed (Kaushansky et al., 1 988; 
Grossman et al., 1989; Demetri and Griffin, 
1 991 ; McKay and Leigh, 1 991 ; Boxman et al., 
1996; Kawada et al., 1997; Mansbridge 
et al., 1999; Mueller and Fusenig, 1999; 
Rennekampff et al., 2000; Gillitzer and 
Goebeler, 2001; Angel and Szabowski, 
2002; Werner and Grose, 2003; Schroeder 
et al., 2009; Yadav et al., 2010). Data are 
means; n = 4 from a single experiment. 
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to ellipsoid shape in basal and suprabasal regions (Fig. 4 D). We 
conclude that basal cells throughout the complete organotypic 
culture collectively change their polarization toward the wound 
edge in all regions of the un wounded tissue in two dimensions: 
in time as well as in distance to the wound margin. Together 
with the fact that the basal cells showed constant adhesion by 
E-cadherin, the collective polarization indeed points to the collec- 
tive migration of the basal cell layer toward the punch wound. 



Basal cells collectively migrate to the tip 
of the EET underneath a mechanically 
stable suprabasal compartment 

The capability of keratinocytes to collectively migrate depends on 
the expression of their junctions. In intact tissue and the EET, we 
characterized the cell junction proteins Dsgl (desmoglein 1), oc- 
cludin, and E- and P-cadherin 48 h after wounding (Fig. 5, A-C 
and F). Generally, we observed cell junctions normalizing quickly 
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Figure 3. Quantification of proliferation by Ki67 image processing of virtual slides. (A) Ki67-stained wound cultures were separated into 1 0 regions of 
identical size equaling five concentric rings. The wound area was divided into four regions (d-g), whereas unwounded areas surrounding the wound were 
divided into three regions (a-c; h-j) at each side of the wound. Broken lines denote the wound margin. (B) The epidermal proliferation index (EPI) and the 
total number of Ki67 + cells were spatially and temporally evaluated by analyzing each region independently. The wound region is highlighted in gray. 
The red dotted lines represent the baseline levels of the epidermal proliferation index. The data shown are from a single representative experiment out of 
two repeats; n = 2 (day 0 [d0]), 3 (day 1 ), 3 (day 2), 3 (day 3), 3 (day 4), 3 (day 5], 2 (day 6), and 3 (day 10]; 113 slides were analyzed in total; data 
are means. Bar, 1 mm. 



toward the intact tissue, whereas the leading edge showed a 
strong P-cadherin expression (Fig. 5 F). All suprabasal cells of 
the EET expressed Dsgl and occludin, indicating intact des- 
mosomes and tight junctions that are immediately formed after 
the EET becomes multilayered. E-cadherin and actin were im- 
munolocalized in basal and suprabasal cells of the EET, show- 
ing functional adherens junctions and microfilaments in the full 
length of the EET (Fig. 5, C and E). The gap junction protein 
connexin 43 was not detected in the leading edge cells but in 
the rest of the EET and intact tissue (Fig. 5 D). Collectively, 
this indicates a mechanically stable suprabasal compartment, 
whereas the basal layer retains its migratory flexibility necessary 
for collective migration (Fig. 5 G). 

Collectively migrating cells continuously 
build a multilayered epithelium in which 
suprabasal cells never contact the ECM, 
thereby ruling out the leap-frog theory 

To understand how the collectively migrating cells form a mul- 
tilayered epithelium, we analyzed keratinocyte differentiation 
(CK10, CK14, involucrin, and filaggrin) in the forming neo- 
epidermis 24, 48, and 72 h after wounding (Fig. 6). 24 h after 
wounding, suprabasal keratinocytes showed a laterally declining 
expression pattern of CK10 and involucrin toward the tip of the 
EET. With progressing time (48 and 72 h), the expression of both 
differentiation markers advanced laterally with extension of the 
tongue. Toward the intact tissue as well as in homeostatic regions 
of the wound models (Fig. 6, H), the differentiation patterns 



were equal to normal human epidermis (Fig. 6, NH). Late differ- 
entiation marker filaggrin was not expressed in the EET between 
0 and 72 h but present in the intact tissue, indicating the time the 
EET needs to further differentiate. At all time points, the EET 
showed a continuously differentiated multilayered epithelium in 
which suprabasal keratinocytes never were in contact with the 
ECM or the basal membrane. 

Cells of the wound margin become 
distributed on top of the EET by 
unstained in-migrating basal cells, 
indicating a lifting mechanism 

We then concentrated on the clarification of how the collectively 
migrating cells form a multilayered epithelium during EET ex- 
tension. To this end, we designed a novel two-step staining using 
the cell trackers green 5-chloromethylfluorescein diacetate 
(CMFDA) and red CMTPX. Although green CMFDA was di- 
rectly applied after wounding, red CMTPX was added 24 h after 
wounding into the punch (Fig. 7 A). This enabled us to track the 
spatial movement of a subpopulation of keratinocytes from the 
wound edge at day 0 (Fig. 7 B, green cells) or the further spatial 
movement of a completely stained EET at day 1 (Fig. 7 B, red 
cells) during wound closure. We showed that both cell trackers 
are not lost over a period of 3 d in primary human keratinocytes 
(Fig. S3 A). Furthermore, the cell trackers had no adverse effect 
on the vitality of keratinocytes in the organotypic wound models, 
proven by MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetra- 
zolium bromide) and TUNEL assays (Fig. S3, B and C). 
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Figure 4. Polarization of collectively migrating basal cells allow shield extension. (A) Histological sections of organotypic wound cultures were separated 
into six different regions dependent on the distance to the wound, and the nuclear major axis of basal keratinocytes was analyzed. The inset shows the 
magnification of the wound margin (box) in region 3 (R3) to illustrate the change of nuclear major axis orientations (red lines). (B) Plotted nuclear orientation 
degrees of the major axis analysis for the respective regions l 2 and 36 h after wounding. (C) E-cadherin staining of representative histological sections 
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Figure 5. Spatial regulation of cell junctions and communication channels in the EET. (A-F) Histological sections of representative EETs showing the expres- 
sion of desmoglein 1 (Dsg 1 ; A], occludin (B), E-cadherin (C), connexin 43 (Cx 43; D), actin (E), and P-cadherin (F). (G) The spatial expression patterns (red 
colored cells) are schematically summarized. Arrowheads indicate the wound margin. Broken lines denote dermal-epidermal junction. Bar, 200 unn. 



At 0 h, the time-lag double stain resulted in the first two to 
three green vertical rows of CMFDA + keratinocytes proximal to 
the wound site showing a homogeneous vertical staining in all 
layers of the section but a very reduced lateral spreading of the 
stain (Fig. 7 B). At 12 h, when keratinocytes had started to mi- 
grate, this cluster of CMFDA + cells was stretched laterally by 
in-migrating CMFDA - basal cells originating from the sur- 
rounding unwounded tissue regions. 

At 24 h after wounding, a long-stretched triangular EET 
had formed with a thin top layer of green CMFDA + keratino- 
cytes. This thin green CMFDA + cell layer connected the granular 
cells of the intact epidermis with the leading edge keratinocytes of 



the EET (Fig. 7 B). In contrast, the basal cells of the EET were 
CMFDA", indicating that unstained cells originating from the 
intact tissue had migrated into the wound region. This in turn 
means that while migrating into the wound, they have lifted the 
CMFDA + cells and thereby formed a multilayered tongue. 

At 48 h after wounding, also, the CMTPX + cells that were 
stained 24 h after wounding now had distributed continuously on 
top of the EET, forming a shield completely covering the wound 
bed. In contrast, the cell layers below this shield were formed 
again by entirely unstained (CMFDA - and CMTPX") cells. 
This showed for the third time (12, 24, and 48 h) how unstained 
basal cells migrated into the wound bed, thereby lifting leading 



of unwounded tissue and wound margin. Compared with unwounded tissue cells within the wound margin show a nuclei movement according to the cell 
boundaries rearward in opposed direction of cell migration. The arrows indicate displacement vectors. Asterisks indicate the geometrical center of the cell. 
(D) Nuclei deformation of wound cultures indicated by white arrows are shown by representative DAPI stains. The arrowhead indicates the wound margin. 
Broken lines denote dermal-epidermal junction. The data shown are from a single representative experiment. For the experiment shown, n = 5; 2,341 cell 
nuclei were analyzed in total; **, P < 0.01, Student's ftest. Bars, 100 unn. 
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Figure 6. Differentiated keratinocytes never contact the ECM while restoring the physiological epidermal differentiation pattern. Histological sections of 
representative EETs, homeostatic wound model regions (H) and normal human skin (NH) 24, 48, and 72 h after wounding. The sections show a double 
fluorescent staining of the differentiation markers cytokeratin 10 (CKIO), involucrin (IVL), and filaggrin (Fil) stained with Alexa Fluor 594 in combination 
with CK14 (Alexa Fluor 488) as a marker for undifferentiated, basal keratinocytes. Arrowheads indicate the wound margin. Broken lines denote 
dermal-epidermal junction. Bar, 300 \im. 



stained cells. To exclude any artifacts caused by loss of staining, 
we performed a quantitative spatial analysis of CMTPX staining 
in EETs at three different time points (0, 12, and 24 h; Fig. S3, 
D-F). We thereby measured the lateral distribution of CMTPX 
and observed gradients with maximum staining intensity at the 
tip of the EET excluding a loss of the dye as a result of prolif- 
eration in the tip. The example of 36 h (Fig. 7 C) further illus- 
trates the rearrangements of cells during wound closure when 
distinguishing three different zones of the EET. The viable 



layers of zone 1 (the wound margin) consist entirely of CMFDA , 
CMTPX" cells. In contrast, CMFDA + , CMTPX + fibroblasts of 
the underlying dermal component serve as a positive staining 
control. In turn, highly differentiated CMFDA + cells in zone 1 in 
close vicinity to CMFDA" cells negate diffusion artifacts. Zone 2 
(the core EET) shows the accumulation of CMFDA + , CMTPX + 
and CMFDA + , CMTPX" cells in the suprabasal and shield layer, 
having previously been lifted by unstained basal cells originating 
from the intact tissue surrounding the wound. Zone 3 (the 
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Figure 7. Capturing the spatial dynamics of the EET by applying a fluorescent time-lag double stain. (A) Experimental design of the fluorescent time-lag 
double stain experiment. CMFDA (green) was subsequently applied, whereas CMTPX (red) was applied 24 h after wounding, and the predicted outcomes 
hypothesized by the leap-frog and sliding theory were compared with the results of the experiment. (B) Sections of the EET showing the initial migration of 
basal CMFDA-labeled keratinocytes 0 and 1 2 h after wounding and the spatiotemporal dynamic behavior of CMFDA, CMTPX-labeled double-positive cells 
24 and 48 h after wounding. The 0- and 24-h images shown here are presented again in Fig. 8 B to represent 1 .5 and 53.3% reepithelialized wounds, 
respectively. The white arrows indicate CMFDA, CMTPX-negative keratinocytes, which have entered the EET from unstained regions. (C) Displays a rep- 
resentative EET 36 h after wounding, which can be separated into three different zones. Zone 1 shows the wound margin where labeled cells had started 
migration. Zone 2 displays the unstained migrating basal keratinocytes originated from tissue adjacent to the wound entering the EET. Zone 3 shows the 
collision of two EETs and the formation of a bulk. Arrowheads indicate the wound margin. Broken lines denote dermal-epidermal junction. The images 
shown are from a single representative experiment out of three repeats; 60 wound-healing models were analyzed in total. Bars, 1 00 urn. 
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collision zone) displays two EETs originating from the left and 
the right wound margin, forming a bulk of stained cells upon col- 
lision, resulting in a thickened epidermis in the collision zone. In 
conclusion, we showed that in epidermal wound healing, collec- 
tive migration delivers cells to the tip of the EET, where they sub- 
sequently become lifted and connected to the extending shield of 
the neo-epidermis, thereby incrementally reconstructing a multi- 
layered epithelium. 

Computational modeling leads to the E5M 

By integration of lifted CMFDA + cells in the EET, an incre- 
mentally extending shield is created (see Fig. 10). We denote 
this the ESM. To test whether the ESM is capable of explain- 
ing the continuous elongation of the EET by collectively 
in-migrating basal cells, we developed a multicellular com- 
putational simulation of wound closure. We used our EPISIM 
platform (Grabe and Neuber, 2005; Sutterlin et al., 2009, 2013) 
to build and simulate the ESM with a computational model 
(see Materials and methods). This model contains four distinct 
cell types. (1) Fast dividing cells: fast dividing cells have a 
fixed position to limit basolateral cell spreading and increased 
proliferative activity to mimic the activated state of the ke- 
ratinocytes under wound-healing conditions. (2) Basal cells: 
these cells were initially marked in green to simulate the cell 
tracker CMFDA and were placed on the basal membrane. (3) 
Suprabasal cells: suprabasal cells are part of the un wounded 
tissue region and function as an initial shield for the migratory 
basal cells. (4) Early suprabasal cells: if a basal cell is lifted 
from the basal into the suprabasal compartment, it differenti- 
ates into an early suprabasal cell. All simulated cells influ- 
ence each other by exerting specific adhesion and intercellular 
pressure forces on neighboring cells (Fig. 8 A and Fig. S4). 
Generally, basal cells exert lower adhesion forces than early 
suprabasal cells or suprabasal cells, thereby mirroring the ex- 
pression of adherens junctions (E- and P-cadherin), leading 
to a flexible migratory stratum basale. The expression of des- 
mosomal cell-cell junctions (Dsgl) in turn characterizes the 
mechanically stable suprabasal cell layers. Impairment of ac- 
tive migration in the wound area is simulated with high basal 
adhesion (see Materials and methods). This reflects within the 
scope of our modeling approach the requirement of leading 
edge keratinocytes to adapt to the new substrate of the ECM 
by an integrin shift as well as the laminin-5 deposition of lead- 
ing edge keratinocytes, resulting in an continuous buildup of 
the basal membrane. 

In the simulation, we modeled cell lifting by (a) the 
pushing force from the proliferative compartment as well as 
(b) the adhesion force between the last shield cell and the 
newly lifted cell (Fig. 8 A). The simulation of our ESM model 
mirrored the cell behavior of migrating keratinocytes in our 
organotypic in vitro wound-healing model (Fig. 8 B, Video 1, 
and Video 2). In the early simulation phase, the EET is visible 
as a monolayer of green CMFDA" 1 " basal keratinocytes mi- 
grating from the wound margin into the wound. This closely 
corresponds to the experimental situation in which the EET 
had covered 8% of the wound area. During simulation, cells 
build up the characteristic triangular structure of the EET, 



showing CMFDA + cells at the tip and in the suprabasal layer 
of the extending tongue. Continuously, unstained basal cells 
migrate into the EET triggering the ESM and further driving the 
wound closure process. The cell distribution obtained at the end 
of the simulation was in accordance to the experimental re- 
sults seen in our organotypic in vitro wound models. The sim- 
ulated tissue kinetics proved the 2D as well as the 3D models 
capability to yield a fully reepithelialized and stratified neo- 
epidermis (Fig. S4 B). 

Leader cells in the 3D situation express 
occludin whose blockage perturbs the ESM 
and results in shedding of bound cells 

Our simulation pointed to two potential causes leading to lift- 
ing of cells: (1) the pushing force from the proliferative com- 
partment as well as (2) cell junctions. Previously, it was shown 
that integrin blocking delays wound closure but does not lead 
to a spatially perturbed EET, indicating that probably cell 
junctions play a more crucial role in cell lifting (Egles et al., 
2010). Our analysis of cell junctions showed that desmosomal 
junctions were not expressed at the tip of the EET, whereas 
adherens junctions were ubiquitously expressed throughout 
the tongue. In contrast to those junctional proteins, only oc- 
cludin showed a colocalization with the cell trackers CMFDA 
and CMTPX (Fig. 9 A). To test whether occludin inhibition 
might be able to perturb the ESM, we blocked tight junction 
assembly during reepithelialization by using a characterized 
synthetic peptide (Wong and Gumbiner, 1997; Nusrat et al., 
2005). This peptide is emulating the occludin loop B, whereas 
its scrambled version is known to still bind occludin, although 
with substantially less affinity. Functional and scrambled pep- 
tides tested in HaCaT cell culture confirmed their accumula- 
tion at the tight junction protein occludin (Fig. S5). 48 h after 
wounding, functional peptides accumulating at the leading 
edge and the monolayer zone of the EET demonstrated that 
leader cells present free occludin binding sites in the 3D situ- 
ation before being included in the shield (Fig. 9 B). At 72 h, 
massive amounts of leader cells had been captured by the 
peptide and subsequently extruded. In contrast, the scrambled 
peptide with low occludin affinity did not completely abolish 
shield extension but led partly to incorporation of peptide- 
tagged leader cells into the extending tongue and to a small 
extent to the extrusion of cells. Unaffected parts of the EET 
as well as the EET of the negative control showed a smooth 
suprabasal layer with no evidence of cell shedding. Thus, 
we conclude that the ESM can be successfully perturbed 
using occludin peptides, pointing to a crucial role of occlu- 
din in wound closure. 

Discussion 

To derive a consistent theory of how the cellular processes of 
migration, proliferation, and differentiation are intertwined in 
skin wound healing, we established a technology pipeline al- 
lowing to analyze 3D tissue reorganization in organotypic in vitro 
wound models comprising a punch wound as well as an intact 
wound margin over a time period of <10 d. 
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Figure 8. Computational modeling of the ESM. (A) The EPISIM platform was used to build a CBM of the ESM. The model contains four distinct cell types: 
(1) fast dividing cells (FDC), (2) basal cells (BC], (3) suprabasal cells (SBC), and (4) early suprabasal cells (ESBC). All cells are placed on a basal mem- 
brane and exert specific adhesion forces, pos., position. (B) During simulation, keratinocytes built up the characteristic triangular structure of the EET, with 
CMFDA + cells becoming distributed on top (dashed arrows), thereby reflecting the experimental results obtained by using the organotypic in vitro wound- 
healing model. Arrowheads indicate the wound margin, and white arrows show the leading edge of the EET. Broken lines denote dermal-epidermal junc- 
tion. The 1 .5 and 53.3% reepithelialized wound images are reproduced from Fig. 7 B (0 and 24 h, respectively). Bar, 100 urn. 



In viva-like organotypic in vitro 
wound-healing model 

The basis of our study is formed by commercially available, 
organotypic skin cultures, developed and manufactured in 
a good manufacturing practice-approved process. In many 
publications regarding different molecular and morphological 
readouts, they have shown a biological response equivalent to 
human skin (Hayden et al., 2009) comprising cytokine secretion 



(Mallampati et al., 2010) as well as proliferative (Black et al., 
2010), inflammatory, and metabolic markers (Hu et al., 2010). 
To further provide evidence of the in vivo relevance of our 
wound model, we measured the temporal profile of key cyto- 
kines regulating cell behavior and compared it with previously 
published clinical skin-wounding data. In human wound fluids, 
the highest cytokines determined were IL-6, IL-8, MCP-1, and 
G-CSF (Grimstad et al., 201 1). In a forensic analysis of human 
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Figure 9. Tight junction perturbation using a synthetic peptide emulating occludin loop B. (A) Histological sections of representative EETs showing occludin 
CMTPX colocalization (white arrows) in the suprabasal compartment over the full length of the EET. (B] A synthetic functional peptide emulating occludin 
loop B was used to block tight junction assembly during tongue extension. The integration of the biotin-tagged peptide into the EET was detected by ap- 
plying streptavidin-FITC (green). Occludin-blocked cells failed to integrate in the ESM and were shed during tongue extension (white arrows). Insets show 
magnified areas of the EET to illustrate the effect of the functional peptide as well as the effect of the negative control. Arrowheads indicate the wound 
margin. Broken lines denote dermal-epidermal junction. Bars: (overview screens) 250 urn; (inlays and insets) 50 urn. 



skin tissue also IL-6 and IL-8 emerged as the most elevated 
cytokines (Takamiya et al., 2008). The same cytokines and 
growth factors have also been reported in wound fluids (Ono 
et al., 1995; Angel and Szabowski, 2002). Our results are in 
accordance to these literature as the proinflammatory cytokines 
IL-6, IL-8, and MCP-1 were found most abundant, showing a 
mean concentration of >20,000 pg/ml. IL-1, which showed a 
peak 1 d after wounding, is one of the key triggers in cytokine 
response (Werner et al., 2007). Our findings therefore point to 
the double paracrine cross talk between keratinocytes and fi- 
broblasts in wound healing (Fig. 2) as proposed by Angel and 
Szabowski (2002). In brief, keratinocytes secrete IL-la, which 
regulates components of the AP-1 transcription factor c-Jun 
and JunB in fibroblasts. This in turn activates AP-l-dependent 
target genes in fibroblasts, leading to a secretion of keratinocyte 
mitogens (e.g., G-CSF and IL-6; Angel and Szabowski, 2002; 
Florin et al., 2004). 



Cell proliferation 

Regarding the dimension of proliferation in the neo-epidermis, 
conflicting information has been reported (Jansson et al., 
1996; Laplante et al., 2001; Patel et al., 2006; Garlick, 2007). 
Our results show proliferative activity of >8 mm around a 
2-mm punch and that the cell supply of the extending shield 
is nearly exclusively generated in the intact unwounded tissue 
and not by tongue extension. Later, during wound closure, 
new cell production is translocated in form of a concentric 
wave toward the center of the punch wound (Fig. 3 D). Hy- 
perproliferation or a proliferative burst at the wound margins 
has been described previously as an initially activating, posi- 
tive feedback after wounding but not as a traveling wave 
(Viziam et al., 1964; Garlick, 2007). Moreover, we precisely 
quantitatively measured this proliferative burst and showed 
its temporal and spatial distribution throughout the EET and 
the intact tissue. 
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Figure 10. The ESM is based on collectively migrating basal cells. (A) Results from the double fluorescent staining experiment show that stained cells always 
become distributed on top of the EET, which can be principally structured in a single layer and a multilayer zone. The broken line denotes dermal-epidermal 
junction. (B) Cells within the suprabasal compartments of the multilayer zone form mechanically stable tight junctions and desmosomes serving as a shield 
for the underlying wound bed. Because of its cell junction composition, this suprabasal layer is more static, whereas the underlying basal compartment 
exclusively expresses E- and P-cadherin, which can rapidly be remodeled, showing a dynamic migration behavior. (C) Because of the mechanically stable 
structure of the shield, extension can only be achieved by adding cells at the point at which the single layer zone connects to the multilayer zone (lifting 
point). At the lifting point, cells are transported upwards by following basal cells, whereas not yet lifted single-layer tip cells continue migration until they 
are lifted by the following cells. (D) Cells of the basal layer that supply the extending shield are in the beginning (day 1 ) produced in the intact unwounded 
tissue by a proliferative burst. With progressing wound closure (day 3), proliferation activity follows the EET to the center of the wound. (E) Viewed from 
top, dynamically moving cell fields can be identified characterized by directed cell migration. Although basal cells from the unwounded tissue show collec- 
tive cell migration to actively supply the shield extension, static suprabasal cells shield the in-migrating basal cells by continuously building a mechanically 
stable barrier. (F) Depending on the contact point of perturbation, it results in a structural disorganization and shedding of occludin-blocked cells or in a 
decrease of reepithelialization speed if integrin p4 is blocked. Bar, 100 urn. 



Collective migration 

The cell junction expression patterns seen in our in vitro wound- 
healing model are in accordance to in vivo wound-healing studies 
(Malminen et al., 2003; Becker et al., 2012; Chavez et al., 2012). 
These patterns indicate a concentric field of basal keratinocytes 
collectively migrating toward the punch wound. During this mi- 
gration, cells remain physically and functionally connected via 
E- and P-cadherin while showing a multicellular polarity (Friedl 
and Gilmour, 2009; Ilina and Friedl, 2009) of the cytoskeleton to- 
ward the punch inside the intact tissue as well as throughout the 
full length of the EET. Polarization was confirmed by computa- 
tional measurement of nuclei rotation and displacement. The ca- 
pability of classical cadherins to control nucleus position and cell 
polarity has been described earlier (Maniotis et al., 1997; Dupin 
et al., 2009). We further included Golgi stainings for final proof 
of cell axis rotation in the intact tissue. Thus, in accordance with 
literature (Maniotis et al, 1997; Friedl et al, 2011), our results 
show basal cells collectively migrating from unwounded regions 
toward the punch wound. The suprabasal compartment thereby 
provides a mechanically robust shield with the described cell 
junctions (see Results). 



ESM for tongue assembly 

So far in literature, three mechanisms have been proposed since 
the last 40 yr that try to explain the recreation of a multilayered 
epithelium: the leap-frog, the tractor-tread, and the Usui model 
(Krawczyk, 1971; Radice, 1980; Paladini et al., 1996; Woodley, 
1996; Usui et al., 2005). Our results show that the EET always 
moves as a triangular multilayered epithelium with full height 
toward the former wound margin and with a single cell tip at the 
protruding end toward the wound bed. Keratinocyte differentia- 
tion markers demonstrated at no time point that differentiated 
keratinocytes got into contact with the ECM, thereby ruling out 
the leap-frog model. We then studied cellular movement of the 
EET by the fluorescent time-lag double stain showing that basal 
cells of a stained EET with progressing time always become 
distributed on top of the EET (Fig. 7, B and C). Our results, 
which are summarized in Fig. 10, lead to the ESM thus ruling 
out all other reepithelialization mechanisms so far proposed in 
literature. Instead, directly after wounding at 12 h, green CMFDA + 
cells are pushed out of the intact wound margin by following 
basally located CMFDA" cells. At no point did we observe 
CMFDA" cells at the tip of the tongue. Thus, the CMFDA + 
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cells formed, in the end, an incrementally extending shield. At 
completion of wound closure, all CMFDA + cells had ended up 
in the outmost cell layer (forming the shield) or still remained at 
the tongue's tip for further migration, as they were not needed 
for shield extension yet. It is thereby remarkable that immedi- 
ately after becoming part of the shield, keratinocytes became 
stretched, in which one single cell shielded frequently more 
than five keratinocytes. By this stretching capability, a rela- 
tively small reservoir of CMFDA + cells from the wound edge 
can shield large neo-epidermal surfaces. 

Multicellular systems biological simulation 
of the ESM 

Our results show that wound healing is a higher level tissue 
function, which only results in an intact tissue again, if the indi- 
vidual cellular behavior is orchestrated in a coherent way. To 
theoretically test the ESM for its capability in providing correct 
reepithelialization and to identify potential regulators, we used 
a multicellular computational simulation. The mechanical pa- 
rameters used in the simulation reflect the adhesion characteris- 
tics of cells in the in vitro wound model, in which basal cells 
showed low and suprabasal cell high adhesion forces. The bio- 
logical motivation for proposing lower adhesion forces in basal 
cells compared with suprabasal cells in our in silico model is 
further supported by the fact that these cell types have different 
mechanical stabilities caused by specific characteristics of their 
cell-cell junctions. In brief, basal cells have small, less orga- 
nized desmosomes with a low electron density (Green and 
Simpson, 2007), which are further destabilized because of the 
low Ca 2+ concentration in the basal layer (Menon, 2002; Proksch 
et al., 2008). The major junction of basal cells are adherens 
junctions binding to actin filaments with a diameter of 7 nm and 
therefore exhibiting less mechanical strength. In addition, basal 
intermediate filaments (CK5/CK14) build tonofilaments, show- 
ing a loose formation with less mechanical stability (Moll et al., 
2008). Furthermore, growth factors such as EGF lead to a dis- 
ruption of hemidesmosomes as a result of a PKC-a-mediated 
phosphorylation of (34 integrins. These integrins dissemble 
from the hemidesmosomes and reassemble on focal contacts, 
thereby leading to an increased motility of basal cells under 
wound healing conditions (Rabinovitz et al., 1999; Litjens 
et al., 2006). 

Suprabasal cells instead have large desmosomes with a 
high organization level and a high electron density (Green and 
Simpson, 2007). These desmosomes have an elevated strength 
caused by a high Ca 2+ concentration in the suprabasal layer and 
are linked to tightly bundled intermediate filaments (CK1/ 
CK10) with diameters of 10 nm, therefore exhibiting maximal 
mechanical strength (Menon, 2002; Niessen, 2007; Moll et al., 
2008; Proksch et al, 2008). 

Our results show that the identified components of the 
wound closure process play consistently together so that, start- 
ing from proliferation, cells are transported toward the wound 
edge and then are assembled in such a way that, according to the 
ESM, a novel 3D epidermis is constructed. Thereby, wound 
closure is achieved with an incrementally growing stable and 
robust shield, under which collective migration supplies novel 



cells for extending the tongue being in turn stepwise displaced 
toward the wound center. 

For realizing the simulation, passive migration through 
displacement by proliferation was sufficient. We explicitly do 
not make further assumptions regarding the actual biomechani- 
cal mode of collective keratinocyte migration in silico as this is 
beyond the scope of this paper. 

The different functions of the apical and 
the basal side of the tongue 

The simulation of the emerging triangular form of the EET 
pointed to the fact that lifting of cells in the ESM toward a mul- 
tilayered epithelium can be regulated either by the lateral or the 
apical side of the tongue. Egles et al. (2010) provided an exam- 
ple of perturbing lateral movement of the EET by blocking of 
integrin (34, which is binding to laminin-5 (laminin-332). Al- 
though their results showed a substantial slowdown in EET 
extension (Egles et al., 2010), it did not lead to the structural 
disturbance of the tongue. To test a perturbation of the apical 
side of an EET, we blocked occludin as an example of a tight 
junction protein by using a peptide emulating occludin loop B, 
resulting in extruded cells during the reepithelialization pro- 
cess. This proves that occludin binding sites are exposed by the 
tip cells and that the free transmembrane occludin binding sites 
are further essential for integrating the in-migrating basal cells 
into the extending shield. Furthermore, this points to a potential 
key role of occludin, as its biomechanical function within the 
EET could be the reason for its duality being involved in cell 
migration as well as in sealing the epidermis during wound 
healing (Niessen, 2007; Du et al., 2010). 

In conclusion 

By revealing the ESM, this work provides a consistent theoretical 
model of how reepithelialization actually occurs in healing skin 
wounds. This shows that wound closure is not a process only hap- 
pening at the leading edge of migrating cells. Instead, closure and 
healing appear as higher-level processes of their own right, on a 
higher abstraction scale than individual cellular processes. 

Materials and methods 

Wound model 

EFT cultures. Human EFT in vitro skin cultures were obtained from MatTek Cor- 
poration. The cultures are derived from human neonatal foreskin tissue and 
consist of normal human epidermal keratinocytes and normal human dermal 
fibroblasts (NHDFs), which have been cultured to form a multilayered, highly 
differentiated model of the human dermis and epidermis. The cultivation was 
performed according to the manufacturer's recommendation. 

Preparation of cellular DEs. For DEs, two different NHDF populations 
originally isolated from >50-yr-old female skin were used. Obtaining skin 
samples was approved by the Heidelberg Ethics Commission. The cells 
were grown in fibroblast culture medium, consisting of DMEM supple- 
mented with 4.5 g/liter glucose, 100 U/ml penicillin, 100 mg/ml strepto- 
mycin, 250 ng/ ml amphotericin B, and 20% fetal bovine serum at 37°C in 
a 5% C0 2 humidified incubator and harvested by trypsinization. NHDFs 
were used within one to three passages. Cells were incorporated into a 
collagen gel prepared by addition to chilled, buffered, pH 7.0, rat tail col- 
lagen I solution (4 mg/ml; Serva). 2.5 ml of the solution was transferred to 
each well of a 6-well plate. Polymerized collagen gels formed after incuba- 
tion at 37°C for 15 min in a 5% CO2 incubator. Once manufactured, the 
DEs were stored, immersed in fibroblast culture medium for 10 d at 37°C 
in a 10% CO2 incubator. During this time, the NHDF secrete additional 
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ECM components for remodeling the collagen matrix and DE, contracting 
thereby ~90%. Medium was changed every 2 d. 

Fabrication of wound cultures. Equilibrated EFT cultures were wounded 
by a circular biopsy punch (Stiefel]. Depending on the respective analysis, 
wounding was achieved either by duplex concentric punching using a 2-mm 
biopsy punch (Stiefel) or by single punching using a 3-mm biopsy punch. To 
quantify reepithelialization and proliferation, the 2-mm punching strategy 
was chosen to generate high amount of data for ensuring in vivo-like wound 
model properties. To analyze keratinocyte migration, a 3-mm punch was 
opted to ensure larger wound space for tracking the newly forming EET. The 
cultures were connected on top to the prepared DE using 6 pi unpolymerized 
rat tail collagen I as "glue" and transferred into 12-well inserts (ThinCert; 
Greiner Bio-One). Cultures were incubated at 37°C for 5 min in 5% CO2 to 
allow polymerization of the collagen glue and terminal bonding of the DE 
with the EFT culture. The manufactured wound cultures were incubated at 
37°C and 5% C0 2 at the air-liquid interface in EFT-maintaining medium 
(MatTek Corporation). Half of the medium was replenished every 2 d. 
Wound cultures were harvested at the indicated time points and subse- 
quently fixed in 4% buffered formalin for 6.5 h at RT. 

Two-step staining procedure 

Keratinocyte migration was observed using the fluorescent cell trackers green 
CMFDA and red CMTPX, which on cleavage of its acetates by cytosolic es- 
terases produce bright fluorescence only in live cells. In total, 60 organotypic 
wound-healing models within three independent experiments were labeled 
stepwise by applying 4 pi of 25-pM CMFDA at time point 0 h or CMTPX at 
time point 24 h directly into the wound and incubation at 37°C for 45 min. 
The wound surface was washed three times with PBS, and fresh medium was 
added in each case. After staining, the wound cultures were incubated as 
long as indicated. 

Histology and automated immunohistochemistry 

Tissue samples were consecutively cut until the wound was microscopically 
visible. For immunohistochemistry, 3-pm sections were deparaffinized by 
reverse ethanol/xylene washes and rehydrated. The histology of the cul- 
tures was assessed by a standard hematoxylin and eosin stain. 

For Immunofluorescent staining, sections were incubated overnight 
at 4°C with primary antibodies, washed 3x with TBS and incubated for 1 h 
at RT with the secondary antibody before mounting in fluorescent mounting 
medium (Dako). The following antibodies were used for immunofluores- 
cence: mouse clone SPM181 anti-filaggrin (1:75), rabbit polyclonal anti- 
laminin 5 (1:1,000), and rabbit polyclonal anti-histone 3 phospho-SlO 
(1:1 ,500; Abeam); rabbit polyclonal anti-connexin 43 (1 :300) and rabbit 
polyclonal anti-P-cadherin (1:200; Sigma-Aldrich); rabbit clone EP700Y 
anti-E-cadherin (1:1 00) and rabbit clone EP1612Y anti-CK14 (1 :50; Epi- 
tomics); clone MIB-1 anti-Ki67 (1:50) and mouse clone HHF35 anti-actin 
(1:1,000; Dako); rabbit polyclonal anti-Dsgl (1:200; Santa Cruz Bio- 
technology, Inc.); rabbit polyclonal anti-occludin (1:75; Invitrogen); clone 
VIK10 anti-CKlO (1:1 00; Acris); mouse clone Sy5 anti-involucrin (1:400; 
Novocastra). Alexa Fluor 488-conjugated goat anti-rabbit IgG and Alexa 
Fluor 594-conjugated donkey anti-mouse IgG (1:100; Invitrogen) were 
used as secondary antibodies. Fluorescent double staining was performed 
simultaneously. For automated immunohistochemistry using the Bond-Max 
system (Leica), sections were incubated for 30 min at RT with primary anti- 
bodies, washed 3x with washing solution provided by the manufacturer, 
and incubated for 8 min at RT with a poly-HRP anti-mouse antibody. The 
following primary antibodies were used for immunohistochemistry: rabbit 
polyclonal anti-GOLGA2 (1:500; Sigma-Aldrich) and clone MIB-1 anti- 
Ki67 (1 :100; Dako). At least two independent consecutive sections of all 
cultures were stained and analyzed. 

Virtual microscopy and image analysis 

Sections were scanned and digitized using a virtual microscopy platform 
(NanoZoomer 2.0-HT; Hamamatsu Photonics). With this system, virtual im- 
ages of whole microscopic tissue sections were obtained. Minor fluctua- 
tions in fluorescence intensity are sometimes apparent in the concatenated 
images but do not reflect differences in image acquisition parameters be- 
tween sections. The imaging systems contains three 4,096 x 64-pixel time 
delays and integration charge-coupled device sensors (cell size 8x8 urn) 
and a 20x objective lens, NA 0.75. Using the NDP.scan software (version 
2.5; Hamamatsu Photonics), slides were scanned at RT in three z layers with 
a spacing of 2 urn each and a resolution of 460 nm/pixel (20x). Images 
were analyzed manually as well as automatically by using the Visiomorph 
image processing software version 3.6.5.0 (Visiopharm). Microscopic 
capturing of whole tissue sections allows for large-scale histological evalua- 
tions with high precision across the complete wound and its outer regions. 



Quantification of reepithelialization. Tissue samples were consecutively 
cut until the wound was microscopically visible. This sectioning plane was 
defined as layer 1 , and 1 0 consecutive sections of 3 urn each were cut. 
Based on that first layer, further defined sections in layer 2 (1 50 pm) and in 
layer 3 (300 pm) were also cut and stored (Fig. 1 D). As equal wound di- 
ameters are a prerequisite for sufficient and comparable reepithelialization 
analysis, we only considered sections having a diameter of 1 ,000-1 ,500 pm 
indicated by the wound margin. Reepithelialization was measured as the 
length of the EET (/eet) in micrometers per day. The maximal wound diam- 
eter (/wound) is indicated by the distance of the unwounded stratum corneum. 
Relative wound closure is obtained by dividing / EET by n/v oun d (Fig- 1 G) As 
the circular punch wound is sectioned, the observed wound diameter /wound 
varies according to cutting depth. EET measurements were performed using 
the NanoZoomer virtual slide viewer software. 

Quantification of epidermal thickness. For quantification of the epider- 
mal thickness, consecutive sections of all three section planes (see Quantifi- 
cation of reepithelialization) within the wound region were used. The basal 
and suprabasal compartment were defined on the basis of the differentia- 
tion markers CK1 4 and CK1 0. For thickness calculation, the area of the re- 
spective compartment was quantified using Visiomorph image processing 
software and divided by the length of the basal lamina. 

Quantification of proliferation (EPI). To analyze the general tissue pro- 
liferation in response to wounding, we quantified consecutive sections of 
all three section planes (see Quantification of reepithelialization) within the 
wound. For detecting and distinguishing negative and DAB-stained Ki67- 
expressing cells, appropriate algorithms were developed and implemented 
using the Visiomorph software. In brief, a preprocessing step using a me- 
dian filter and a Bayesian classifier were applied for nuclei detection. The 
following postprocessing steps have been performed within the hematoxy- 
lin channel to identify nuclei: (a) Small particles beyond 3-pixel particle size 
have been removed, (b) Structures showing a circularity value of <2 were 
identified as potential nuclei, (c) Structures having an area >60 pixels were 
marked as potentially conglomerated nuclei and treated with an eroding 
step of 2 pixels to separate these objects, (iv) An additional separate-object 
step was applied to insufficiently separated conglomerated nuclei, (v) Sepa- 
rated conglomerates were revised and again tested for size and circularity, 
thereby removing objects showing a circularity >2.5 or a size <4 pixels. 
Verified nuclei were identified as proliferative if the Ki67-indicating DAB- 
staining was in a 1-pixel distance to hematoxylin-indicated nuclei. To calcu- 
late the EPI, the amount of Ki67-expressing basal keratinocytes was divided 
by the whole number of basal keratinocytes, to determine the percentage of 
proliferating cells as an indicator for proliferative activity. Furthermore, the 
total number of Ki67 + cells was determined. For this purpose, we harvested 
46 tissue samples and determined Ki67 by analyzing 1 1 3 tissue sections in 
total. For spatial analysis of the tissue sections, virtual slides were separated 
into 10 different regions as follows: The wound area was separated into 
four regions (Fig. 3 A, d-g) of equal size, whereas the unwounded tissue 
surrounding the wound was separated into six regions of equal size, three 
regions (Fig. 3 A, a-c and h-j) of each side of the wound area (Fig. 3 A). 
The EPI and the total number of Ki67 + cells were calculated independently 
for each region (day 1-7). 

Cell polarization 

To assess nuclear shape, we separated the cultures into six equal regions 
(R1-R6) and analyzed every region independently. The epithelial and 
basal regions were manually marked using roipoly function in MATLAB 
(MathWorks, Inc.; for source code see Supplemental material), and the nu- 
clei in the respective regions were detected in the DAPI channel. After ini- 
tial cleaning and intensity adjustments, the nuclei were detected by stepwise 
thresholding. The detected nuclei clusters were further filtered by size and 
split using a watershed method with local maxima as seed points. Ringlike 
stained nuclei were excluded from the watershed procedure. Finally, valid 
nuclei were selected based on region properties such as size and solidity. 
To show the tilt in the nuclei with respect to the basal lamina, the angle of 
major axis was calculated. The orientation angle of a nucleus is measure 
of the angle (in degrees ranging from 90 to 0°) between the x axis (basal 
lamina) and the major axis of the ellipse that has the same second moments 
as the nucleus (Fig. 4 A, red line in nuclei scheme). The regions of eight 
samples were used for the calculations, and in total, 2,340 basal cells 
were analyzed. 

Multiplex analysis of cytokine concentrations 

Multiplex analysis of cytokine concentrations were conducted using the Lu- 
minex 1 00 platform built on xMAP technology (Luminex Corp.). The princi- 
ple of this assay is similar to a capture sandwich immunoassay: By coupling 
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an antibody directed against the target protein to internally color-coded micro- 
spheres, the target protein is bound to the individual microsphere and detected 
by a secondary biotinylated antibody. Streptavidin-phycoerythrin is then used 
to detect secondary antibodies. Based on this principle, the technology en- 
ables measuring of multiple total target proteins within one well. Supernatants 
of wounded and unwounded organotypic cultures before (day 0) as well as 
1 and 3 d after wounding were collected, and concentrations of the following 
cytokines were analyzed: IL-1 a, IL-3, IL-6, IL-8, MCP-1 , and G-CSF. 




Accordingly, the optimal distance between two ellipsoidal cells is 



Computational modeling using the EPISIM platform 

Graphical cell behavioral models (CBMs) are built with the EPISIM graphi- 
cal modeling system. The graphical modeling system allows formulating 
CBMs as graphical process diagrams. Such a process diagram can include 
deterministic as well as stochastic model elements (Sutterlin et al., 2013). 
Thus, pure deterministic, pure stochastic, and also hybrid models can be re- 
alized. A graphical CBM is automatically translated in highly efficient exe- 
cutable code by the EPISIM code generator. This code is then loaded by the 
EPSIM simulation environment, conducting a multiagent-based multicellular 
tissue simulation based on the translated executable CBM (Sutterlin et al., 
2009). The computational model of the ESM, the EPISIM modeler, and the 
EPISIM simulator can be downloaded by using the following link: http:// 
www.tiga.uni-hd.de/downloads.html. 

Cell center-based biomechanical model. Beside the CBM, such a tis- 
sue simulation is composed of a biomechanical model (BM). The BM is 
an off-lattice, center-based model, allowing cells to move freely in space. 
Cells are modeled as discrete objects, and all occurring forces act on a 
cell's center of mass (center based; Dallon, 2007; Drasdo, 2007). In such 
a model, cells are commonly represented by circles and spheres (Walker 
et al., 2004; Galle et al., 2005; Schaller and Meyer-Hermann, 2007; 
Adra et al., 2010) or ellipses and ellipsoids (Palsson and Othmer, 2000; 
Dallon and Othmer, 2004). We developed a 2D as well as a 3D cell 
center-based BM with an elliptical cell shape in 2D and an ellipsoidal cell 
shape in 3D. As depicted in Fig. S4 A, the model aims at optimizing the 
distance between a cell c and its neighboring cells c n ; - This optimization 
is performed for all cells in the multicellular tissue simulation and results in 
passive migration of the cells. A cell's position is equivalent to its center of 
mass, r. The basal membrane is discretized into sections of equal length. 
The center of such a section is denoted by r bm i - Fig. S4 A shows the op- 
timal distance balancing between two consecutive simulation steps t — 1 
and f, considering cell c colored in red. The distance balancing involves the 
calculation of intercellular pressure (see cell c ni ) as well as cell-cell adhe- 
sion (see cells c n 3 and c n 4 ) and finally adhesion to the basal membrane. 
The calculated balance of cell distances is perturbed by cell division. A 
proliferating cell grows and compresses its neighboring cells. These cells in 
turn passively migrate away from the dividing cell because of this intercel- 
lular pressure. Further parameters of the BM are listed in Fig. 8 A. 

Optimal distance calculation. The optimal distance between two el- 
liptical or ellipsoidal cells is the target value of the cell distance balancing 
process, calculated over multiple simulation steps. An ellipse's semimajor 
axis is denoted by a, and the semiminor axis is denoted by b. Accord- 
ingly, the semiprincipal axes of the ellipsoid are a, b, and c. An ellipse's 
or ellipsoid's axes are aligned with the axes of the used Cartesian coor- 
dinate system. 

The calculation of the optimal distance between two ellipses as 
well as two ellipsoids follows the same scheme: (a) calculation of the line 
connecting cell centers (of mass) r c and r„ as well as the same line in the 
opposite direction (r„ and r c ), (b) calculation of the intersection points be- 
tween these lines with the ellipses or ellipsoids, and (c) summation of the 
length of the line segments between a cell's center and the according in- 
tersection point. Consequently, the optimal distance between two cells 
represented by ellipses located at r c and r„ is 
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Intercellular pressure and cell-cell adhesion. The optimal distance is 
used to calculate the intercellular pressure F pr (r c , r„) between overlap- 
ping cells: 



fpr( r c' r n) : 



[ 8 ol xcf opt( r c' , n) . f 



with v nc = r c — r n . The factor 8 0 | (0 < S Q | < 1 ) modulates the allowed 
overlap between two adjacent cells and, by that, a cell's compressibil- 
ity. We used 8 0 | = 0.95 corresponding to an allowed cell overlap 
of 5%. 

The adhesive force F a dh(r c , r„) between adjacent cells is 



F adh( r c r n) = 



l c/ opt( r c- r n) 
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if r„ * r n 
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(1) 



Thereby, cell-cell adhesion is limited to cells whose distance is below 
8 a dh x c/ op i(r c , r„), with 8 a jh = 1 .25. Intercellular pressure F pr (r c , r„) and cell-cell 
adhesion F 0 dh(r c , r n ) are used to calculate a cell movement vector: 



(2) 



- ad \„(l-Fadh(r c ,r n ))xv, 
0 else 



f F pr (r c ,r n )>l 

f F pr (r c ,r n )<lAF adh (r c ,r n )>0. 



The adhesion coefficient adh /c cn can be interpreted as a linear spring constant 
modulating the strength of the adhesive force between two cells. The adhe- 
sion coefficient can be used to model different cell-cell junctions, such as in- 
tegrins, desmosomes, or tight junctions. 

For each adjacent cell, a movement vector m(r c , r„) is calculated. All 
vectors are finally summed and lead to a cell's final position. 

Adhesion to basal membrane. The basal membrane is discretized into 
segments of equal length (10 urn per segment) to model basal adhesion. 
A segment's center is denoted by rb m . The adhesive force F a di,(r c , r„) and the 
cell movement vector m(r c , rb m ) are calculated according to Eqs. 1 and 2. 
The used optimal distance between a cell and the basal membrane is for 
elliptical cells 
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and for ellipsoidal cells 
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The strength of the cell junctions to the basal membrane and the adhe- 
sion, respectively, depends on the adhesion coefficient adh k cn and the cell 
contact time f con _b m . The cell contact time corresponds to the period of time 
a particular basal membrane segment (BMS) was covered by a cell. The 
higher the cell contact time, the darker the coloring of a segment (Fig. S4 A). 
Moreover, the lower the contact time, the less ECM is assumed to exist on 
this BMS. The absence of ECM, in turn, hinders active cell migration. It is 
a property of the chosen modeling approach that cells can migrate only 
passively. Because of this limitation to passive migration, the inability of 
active cell migration is mimicked by high adhesion on the BMSs with low 
ECM and accordingly low f C on_bm values. This high adhesion biologically 
corresponds to low adhesion and to the force feedback of cells that can- 
not actively migrate. We therefore define a high basal membrane adhe- 



sion factor odl, /cc_bm_hi- The effective adhesion coefficient a ^ L - 
the ECM level on a particular BMS is 
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Threshold T con bm is the contact time assumed for reaching the maximum 
ECM level on a particular BMS. The range of the effective adhesion coeffi- 



, according to Eq. 3 is k c _b m 
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Peptide perturbation 

For tight junction perturbation, a functional and a scrambled peptide were 
used showing the following amino acid sequences: Biotin-SQIYALCNQ- 
F YTPAATGLYVD-N H 2 (functional peptide) and Biotin-TSPYAQIYIANFQD- 
TALGYCV-NH2 (scrambled peptide). For 2D experiments, HaCaTs were 
cultured using culture slides until they reached a confluency of 70-80%. 
Cells were washed with HBSS^~, and tight junctions were disrupted by 
using a Ca 2+ switch protocol according to Nusrat et al. (2005). Peptides 
having a concentration of 200 pM were applied, and cells were incubated 
at 37°C for 5 and 24 h. Organotypic cultures were wounded, and 4 ul 
peptide solution having a concentration of 400 uM was applied every 
hour directly into the wound for a period of 6 h. 

Online supplemental material 

Fig. SI shows the comparison of the EPI to the mitotic cell fraction and the 
morphometric analysis of the epidermal thickness during reepithelializa- 
tion. Fig. S2 shows the Golgi positioning in polarized migrating cells. Fig. 
S3 shows the persistence of the cell trackers CMFDA and CMTPX in primary 
human keratinocytes and their influence on the viability of organotypic cultures. 
Fig. S4 shows the in silico biomechanical cell model of the ESM. Fig. S5 shows 
the characterization of functional and scrambled peptides emulating occludin 
loop B on primary keratinocytes and organotypic cultures. Video 1 shows the 
2D simulation of the ESM. Video 2 shows the 3D simulation of the ESM. The 
MATLAB source code of the major axis analysis is provided online as a Word 
document. Online supplemental material is available at http://www.jcb.org/ 
cgi/content/full/jcb.201 2 1 2020/DC1 . Additional data are available in the 
JCB DataViewer at http://dx.doi.org/ 10.1 083/jcb.20 1 2 1 2020.dv. 
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